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Using the Bethe ansatz analysis as was reformulated by Edwards, we calculate 
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the spectral properties of a particle interacting with a bath of fermions in one di- 

m . 

■ mension for the case of equal particle- fermion masses. These are directly related to 



singularities apparent in optical experiments in one dimensional systems. The or- 



■ thogonality catastrophe for the case of an infinite particle mass survives in the limit 

of equal masses. We find that the exponent (3 of the quasiparticle weight, Z ~ N $ 
is different for the two cases, and proportional to their respective phaseshifts at the 
Fermi surface; we present a simple physical argument for this difference. We also 
show that these exponents describe the low energy behavior of the spectral function, 
for repulsive as well as attractive interaction. 
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I. INTRODUCTION 



The motion of a particle interacting with a Fermi bath has been the subject of many 
studies as it relates to a variety of physical phenomena : muon diffusion in metals |l| . 
stability of the fully ferromagnetic state (the Nagaoka problem) 0] and particularly relevant 
to optical experiments in degenerate semiconductors and metals |3,|J], where the hole created 
by optical excitation plays the role of the particle interacting with conduction electrons. 

In the late 60's the Fermi edge singularities in the X-ray absorption spectra of metals were 
explained in terms of correlation effects between conduction electrons and the photon created 
hole; Mahan first predicted the singularities before Nozieres and DeDominicis calculated the 
exact spectrum; these pioneering works gave rise to the so called MND problem ||. Relying 
on perturbative arguments valid for a three dimensional system, it was believed that such 
singularities would disappear in valence band excitations where the mass of the hole is 
comparable to the conduction electron mass, in contrast to the case of core level excitations, 
where the mass of the localized hole can be considered infinite 0. In lower dimensions the 
situation is different || and in fact strong enhancement at threshold was recently observed 
in photoluminescence on quantum wires These experiments have raised further interest 
in optical singularities in low dimensional systems. 

This model was also studied in connection with ferromagnetism in single band models 
and it has recently gained renewed attention in two dimensions because of the possible 
relevance of the single band Hubbard model to describe the high-Tc superconductors. In 
this context it describes the dynamics of a single spin flip in a ferromagnetic background and 
adresses the question of the stability of the Nagaoka state or the limitations of Gutzwiller 
type variational wave functions used for this problem 0. 

In this work we study the spectral properties of the particle in one dimension and for the 
special case where the mass of the particle is the same as that of the fermions in the bath. 
For this analysis we use the analogy to McGuire's solution of the problem of a single spin 
flip interacting with a ferromagnetic background in a one dimensional continuous system, 
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considering the particle as the reversed spin |7]]. Actually McGuire's work is a precursor, 
for the continuum model, of Lieb and Wu's Bethe Ansatz (BA) exact solution of the Hub- 
bard model in one dimension ||. McGuire calculated static correlation functions as well 
as the effective mass of the reversed spin; a lattice version of the latter was presented in 
Ref. ||. However, despite the BA solution, an exact calculation of the dynamic quantities 
has not been possible as it involves the calculation of matrix elements between, difficult to 
handle, BA wavefunctions for the excited states. In this work we use a new presentation of 
McGuire's solution, due to Edwards ||10|| ; the relative simplicity of the wave functions in this 
case allows us to evaluate directly the spectral weight Z and the spectral function A(k, to) 
in the small momentum k regime. Our results show that the orthogonality catastrophe of 



Anderson occurs [|Tl|] : a quasiparticle description of the excitations in terms of non inter- 
acting eigenstates is therefore no more possible; the spectral function is totally incoherent 
and has the same sort of divergence as in the infinite mass case. 

The occurence of the orthogonality catastrophe and the divergence of the spectral weight 
were recently predicted within perturbation theory for the one dimensional problem ||. The 
edge singularities of the absorption spectrum were also derived for a Tomonaga-Luttinger 
model and it was claimed that the corresponding exponents would not depend on the mass of 



the hole ||12|| . In our work however, we show that the exponents of the spectral function differ 
in the two extreme limits of infinite mass and equal masses. This seems in contradiction 
with the latter results; we will discuss the importance of backscattering, which was omitted 
in the Tomonaga-Luttinger model, to correctly describe the case of an infinite particle mass. 
We note that, except for the very special case of one hole in a half filled band with infinite 



repulsion ||13|| , our work presents the first calculation of dynamical correlation functions using 
directly the BA wavefunctions. Our calculations agree with previous results by Frahm and 
Korepin |14[ thus giving support to the assumption of conformal invariance of the model 



used. 

In summary we show in this article that the orthogonality catastrophe takes place for the 
case of equal particle-fermions masses and that the same exponents characterize the diver- 



gence of the spectral function at threshold in the small momentum regime. The exponents 
however differ in the two cases of infinite mass and equal masses. The article is organized as 
follows. In section 2 we present the model and its exact solution following mainly Edwards. 
Our analysis is performed for the continuous model and we simply quote the results for the 
lattice. In section 3, we derive the orthogonality catastrophe for the repulsive interaction. 
In section 4, we calculate the spectral function for the repulsive and attractive case; in the 
latter we discuss the influence of the bound state. 



II. MODEL 

The model Hamiltonian describes N fermions of mass m and one particle of mass rrih 
interacting via a delta function potential and moving on a line of length L with periodic 
boundary conditions : 

1 N d 2 Id 2 N 

Throughout this article we use the convention m = 1 and rrih — > rrih/m. When rrih — > 
oo the problem reduces to a single particle problem. The fermions evolve in the static 
potential created by the hole. The response of the Fermi sea to the sudden appearance of 
the disturbance was calculated by Nozieres and De Dominicis [||] . For a finite mass however, 
this model gives rise to a full many body problem. 

We will present at first the problem quite generally for an arbitrary mass and then focus 
on the solution for rrih = 1 following mostly Edwards presentation. It is convenient to 
express the problem in the reference frame attached to the particle. The wave functions in 
the two reference frames are related by a simple translation : 

*(ar , ...,x N ) = e iQxo f(x! -x ,...,x N - x ) (2) 

Q is the total momentum. For our periodic boundary conditions Q = 2iim/L, m being an 
integer, and jjj = Xj — x G [0, L]. The Schrodinger equation for f is then : 
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I 1 / N f) \ 2 i N a 2 N 

It describes the motion of A" fermions in a static local potential and interacting via their 
total momentum. In order to satisfy the equation everywhere and the periodicity, we add 
the boundary conditions : 



L - o 

Vi=o - u 



f(Vi, ■ ■ ■ ,Vn) 
1 d 

— j r f(yi,---,yN) J. =0 = 2Uf{ yi ,..., yi = 0,...,y N ) (4) 

where \x = + mj,) is the relative mass. Edwards proposed the following solution in 

the case = 1 : 

f{y u ...,y N ) = -^=L= det(0 3 -(y f )) (5) 

For this wave function, the cj)j are normalized and the boundary conditions imply : 0^(0) = 
4>j{L) and 0^(0)— 4>j(L) = U<fij(0). This can be achieved if the functions 4>j have the following 
form : 

N 

Hy) = H^ lkty (6) 

t=0 

where the momenta kj are all different and solutions of the BA type of equations : 
Lkj = 2ir (rij + 1/2) - 2arctan(2(A; i - A)/U), j = 0, N 

N i TV 

3 = ^- and i? = -E^ (7) 

j=0 A j=0 

These are Lieb and Wu's equations in our special spin sector of a single spin flip from the 
totally ferromagnetic state and in the continuum limit; for an even total number of particles 
N + 1, the rij are integers. A similar set of equations was originally derived by McGuire. 
The eigenstates are specified by the quantum numbers rij, j = 0, . . . , N and m; A is a real 
number which ensures that the total momentum is Q = 27rm/L for a given choice of rij. 

We would like to focus on the simplicity of this solution compared to McGuire's original 
one. Indeed the latter was written in the static reference frame and defined by pieces in 



the different regions Xj < xo < xy corresponding to a given ordering of the particles on the 
line; moreover it was a superposition of N + 1 determinants. In contrast Edwards solution 
is defined on the compact domain [0, L] N and it takes the form of a single determinant. 
However, we point out that the plane waves appearing in (6) are not the usual free states 
(as k t 7^ 2-kh/L for U ^ 0) and are not periodic on [0, L] unless U = 0. But the functions 
(f)j have to be periodic. 

In order to construct the total wavefunction, we have to first solve the BA equations for 
a given choice of quantum numbers, form the linear combinations in (6) and then build up 
the determinant. We have N + l plane waves and N (fi/s. Is the total wavefunction uniquely 
defined by this procedure ? In fact there is an additional constraint on the coefficients a* 
coming from the periodicity requirement : 

£>$ (l-e tktL ) =0 (8) 
t=o 

This restricts the (fi/s to a N dimensional subspace where we can make an arbitrary choice 
of basis. The total wavefunction being a determinant is thus unique. 
We can find a physically convenient basis of functions : 



.i 



(y) = - sin(^) <?(</)), j = 1, . . . , N 

g(y) = ^ . m (9) 

where the phaseshifts are given by the BA equations : 

n (2k — A\ 

Lkj = 2-KTLj — 25 j, 5j = — — + arctan ( — 3 — — j (10) 

g(y) is an almost localized function for large L and iV; it corrects the plane wave around 
the origin in order to achieve periodicity without affecting its plane wave character almost 
everywhere. We show in the appendix the following statements : 

4-->i(l + 0(l/L)) (11) 

J%]{y)Hy)dy = h + 0{i/L) (12) 



Xf = A= f L e-*^i(y)dy = ( 8m [ 5j \ , + OQog(L)/L), p integer (13) 
J VL Jo [rij -p) - bj/n 

These last results tell us that the <pj behave like an orthonormal set of single particle scat- 
tering states for large systems : 

h(y) ~ l=S k >v + V , 0«y«L (14) 
VL 

These are central results in our paper and the overlaps Xj are used in all the subsequent 
calculations. 



III. ORTHOGONALITY CATASTROPHE 



Anderson studied the influence of a static potential on a Fermi sea [Oj]; he showed that 
the ground state, |^o), of the system in the presence of the potential became orthogonal in 
the thermodynamic limit to the ground state without potential, \ipo) : 

Z H<Wo)| 2 ociV-^-, ^=(f) 2 (15) 

where bp are the phaseshifts of the electrons at the Fermi surface for the two channels 
of even and odd parity wave functions, N the number of fermions and the overlap Z the 
spectral weight. This orthogonality catastrophe takes place in the MND problem where the 
static potential is created by the core hole. Moreover the same exponents f3 ± appear in the 
spectral function of the hole. In our model with — > oo, the odd states do not feel a 
potential located at the origin and the odd phaseshift vanishes; only the even parity states 
contribute to Anderson's orthogonality catastrophe. 

In this section we calculate the spectral weight for our model with equal masses and a 
repulsive interaction. We show that the orthogonality catastrophe also takes place but with 
different exponents than for an infinite mass. The spectral weight can be evaluated with the 
use of our wave function for the ground state and the corresponding overlaps X™ : 

dx ... dx N f* f) 2 = det(A7™) 2 (16) 



/ is here the interacting ground state wavefunction in the reference frame of the particle 
and / the corresponding wavefunction for the non interacting system. The distribution of 
momenta for the ground state is given by the following choice of quantum numbers : the rij 
are consecutive integers from — (N + l)/2 to +(iV — l)/2 and A = 0; the total momentum 
vanishes. Both k and —k appear in the spectrum and it is natural to form even and odd 
combinations of the interacting single particle wavefunctions fy; they are then combinations 
of cos (kjy + 5j) and sin (kjy + Sj) for the positive kj. We note here that both the even and 
odd channels have a non vanishing phaseshift in contrast to the infinite mass case. The 
non interacting states can also be expressed as odd and even functions. The spectral weight 
factorizes then into two contributions coming from both parities, Z + and Z~ : 



The indices I, n, n' run over integers from to (N — l)/2 corresponding to the positive 
solutions of the BA equations. Anderson calculated this determinant using an algorithm 
due to Cauchy; his result reads : 



Sp is the phaseshift at the Fermi surface which can be expressed in terms of the density of 
states n F = l/nk F : 



The exponent for the lattice has the same form but with the corresponding density of states. 
We can now compare our result with the known result for the static impurity. So in the 

rrih = 1 case, we have : 




(17) 



Z ± oc N- (5f/7t)2 



(18) 




(19) 



Z oc N'^-P where 




(20) 



In contrast in the infinite mass case : 




(21) 
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For mh — > oo, the only contribution to the exponent comes from the even parity, the odd 
phaseshift vanishing. For = 1, both parities contribute equally; moreover, the density 
of possible excitations at the Fermi surface is reduced by a factor of 2. This difference can 
be understood easily by simple perturbative arguments. For an infinite mass particle, both 
forward and backward scattering at the Fermi surface involve low energy excitations; for 
m/j = 1 however, backward scattering is forbidden because the recoil energy of the particle 
is of the order of the Fermi energy. So only half of the excitations are allowed. In figure 1 
we show the exponent for the two cases of equal masses and infinite mass as a function of 
the dimensionless parameter Unp. 



IV. SPECTRAL FUNCTION 

In this section we use our BA eigenstates to evaluate the low energy part of the spectral 
function A(Q,u); it is defined as the Fourier Transform of the propagator of the particle 
G Q (t) : 

1 f°° 
A(Q,u) = -Im e wt G Q {t) 

7T J-oo 

G Q (t) = i < T {b Q (tp Q (0)} > (22) 

bq ( resp. bq) is the creation (annihilation) operator of the particle in a free state of 
momentum Q. The propagator can be expanded in the eigenstates of the interacting system : 

G Q (t) =T,^ {Ex ^ o)t \(MQ)\MQ))\ 2 o(t) (23) 

A 

\ipo{Q)) is the non interacting ground state of N fermions with the particle at momentum 
Q; its energy is E Q . \ip\(Q)) are the eigenstates of the interacting system with energy E\ 
and total momentum Q. 

The spectral function describes the photoemission spectrum in the so called intrinsic 
approximation For an infinite mass particle it has no momentum dependence and is 
totally incoherent |TB[ : it has no quasiparticle peak because of the vanishing spectral weight 
and presents a power law singularity : 

9 



In this section we show that the spectral function for the = 1 has the same low energy 
behaviour for momenta close to the bottom of the band Q ~ 0; we treat the case Q = 
in detail and mention the results for other total momenta. We examine first how the low 
energy spectrum of the interacting system with Q = can be described in the language of 
the non-interacting system. 

We are going now to discuss separately the case of a repulsive (U > 0) and attrac- 
tive interaction (U < 0); the difference is the presence of a bound state for the attractive 
potential. 

A. U > 

There are two types of excitations : 'single particle' excitations where one of the quantum 
numbers rij is changed and collective excitations where A is changed. In general, we must 
combine both types to have an eigenstate of zero total momentum. We begin our discussion 
with the simplest excited states that are exactly single particle like and then extend our 
description to the others. 

For the states with a symmetric distribution of the half integers rij + 1/2 around the 
origin, we can achieve a zero total momentum with A = and the distribution of momenta 
is symmetric as well. Moreover when A = 0, the BA equations are decoupled and the 
excitations are strictly single particle like, in the sense that the spectrum corresponds to the 
ground state momentum distribution plus a symmetric momentum distribution of 'particle- 
hole' like excitations. 

The situation is more complicated for the states with an asymmetric distribution of 
quantum numbers because zero total momentum implies A ^ 0. This causes a global shift 
of all the momenta compared to the ground state and these excitations are not strictly 
single particle like. However the excitation energy and spectral weight can be approximated 
in order to recover the single particle description. 
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Let's take the simplest type of these A ^ excitations, namely when one of the rij close 
to the Fermi surface is changed : n'- — rij = J. If J is small compared to N the excitation 
energy is small and A is small as well. A similar analysis of the spectrum as in Ref. J/] can 
be performed; the main results are as follows : 



m - h) 



i 
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A= 



Ul + (2k/Uy U 3 l + (2k t /U) 2 
2vrJ 7T 



1*3 



(25) 



L 2arctan(2A; F /[f) 

In order to evaluate the propagator, we need to know the excitation energy and the 
spectral weight. The change in energy is : 



E\ — En 



2vrJ 



2vrJ 



m* \L Z 



(arctan(2fc F /f/)) 2 



(26) 



The excitation energy E\ — Eq is a sum of two contributions, the first corresponding to 
the particle-hole excitation and the second to a rearrangement of the Fermi sea, a backflow 
term. As pointed out by McGuire the momentum distribution of the particle is centered 
around K = 2nJ/L so that we can interpret the usual backflow as a recoil of the particle 
to the particle-hole excitation in the Fermi sea; its mass is renormalized to m*. For low 
lying excitations E\ — Eq -C l/2k\ so that \2ttJ/L\ kp. In that case the recoil energy is 
negligible and the excitation energy is of particle-hole type only. 

The shift of all the momenta (ki —>■ k[) influences also the spectral weight; as for the 
ground state, it is a simple determinant of the individual overlaps Xf . Because of the 
particle-hole excitation rij — > n'j, the corresponding overlap Aj is replaced by X? ; this is 
the usual effect of an excitation in an independent particle problem. In addition, all the 
other overlaps are influenced by the backflow. Using (25) we find for j ^ I : 



Xf - Xf 



2%\J\ 



Lk r 



1 



2(n,-p)-l 



pi 



xf 



2tt|J| 
Lkp 



< 1 for ki/U -> and k t ~ k F 



< 1 for ki/U -> oo 



(27) 
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Thus the spectral weight has also the usual independent particle form : if one particle-hole 
pair is created around the Fermi surface, only one of the overlaps is significantly altered. All 
the low lying excitations are additive in the sense that their excitation energy and momentum 
is a sum of individual contributions. This allows us to generalize our results to all low lying 
excitations. 

In summary we note that the backflow is not important in our problem, neither for the 
excitation energy nor for the spectral weight. This is simply due to the fact that for a 
scattering between the particle at the bottom of the band and one electron at the Fermi 
surface, for small momentum transfer q, the recoil energy E = q 2 /2 is negligible compared 
to the particle hole energy which is linear in q. We stress that this is no more the case when 
the particle does not lie at the bottom of the band. 

Now the propagator can be approximated in the following way : 



2 



G (t) ~ J2 e^( fe HH(*+*0/2 ( det 

{ni} 

= det f f; e ^-P P ){t+i^)/2 X v x A (28) 
y=-oo / 

where we introduced an energy cuttoff £o corresponding to the range of validity of our single 

particle description. The momenta k p are the non interacting one and the ki are solutions of 

the BA equations for A = ie when we neglect the backflow term. In this last expression we 

can again use a basis of definite parity; this was not possible for the exact eigenstates but 

our description of the low lying excitations allows us to recover this symmetry like in the 

ground state. The propagator factorizes then in two equal contributions for odd and even 

parity. The calculation of this determinant was performed by Nozieres and Combescot and 

we only quote here their result JT7|]: 



GoW-ie-fxoiiZotr^-f-Qit) (29) 

where f3 + = f3~ = (S(cf) / tt) and 5 is a determination of the phaseshift at the Fermi surface. 
Therefore the spectral function has a divergence at an energy loq = —(2/tt) Jq f S(e)de and 
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the exponent is 1 — j3 + — f3~ . The determination 5 is a priori unknown, but if we compare 
the threshold energy ujq with the correct energy shift due to the interaction E$ — E$ = 
— (2/7r) Jq f 5(e)de we can fix the determination : 5 = 5. In summary, the spectral function 
at the bottom of the band has a power law singularity and we recover the exponents f3 + and 
P~ of the orthogonality catastrophe. 

What changes if the momentum is not strictly zero but still much smaller than hp ? The 
energies are simply shifted by a recoil term of the particle Q 2 /2m*; the overlaps Xf around 
the Fermi surface are not affected as well and the form (29) of the propagator remains : 

G Q {t)~e- lQ2t / 2m *G {t) (30) 

The spectral function is then only shifted rigidly by the recoil energy : its exponent is 
the same. Nevertheless this is true only close to the bottom of the band : the numerical 
calculations showed that this shape was strongly altered in the vicinity of the Fermi surface 

0- 



B.U<0 

What changes for the attractive potential ? The main difference comes from the appear- 
ance of a bound state. Indeed the BA equations have two imaginary solutions which in the 
large L limit are : 

K = A±i— (31) 

They always appear in pair in order to have a real total momentum. The ground state is 
reached when A = 0; the other momenta are solutions of the BA equations with U < 
and the integers rij are consecutive from — (N — l)/2 to + (N — 3)/2. We would like to 
note that the phaseshifts do not behave like in 3D : the appearance of a bound state is not 
characterized by a phaseshift going to tt at the bottom of the band but it is simply opposite 
to the phaseshift for the repulsive potential with comparable strength \U\. 
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The ground state wave function is simple as well : it looks similar to the repulsive case but 
with one of the (fi/s describing a bound state : 

Mx) = A cosh (^(x - |)) , A oc e- UL ' 2 (32) 

All the analysis is similar to the repulsive case; we should simply add one bound state 
in the single particle like excitation; its overlap with the free states is simply : 

= z QB (33) 

B stands for bound state, e# = k 2 /2, «b = cos\i(kL / 2) k, / \^A~o and e m is the free state 
energy. 

Nozieres and Combescot studied the influence of the bound state and they showed that 
the spectral function had two divergences, one at the energy of the true ground state, the 
other at the energy of the lowest excited state which does not contain any bound state. We 
first consider the absolute threshold. 

We can perform the same analysis of the low lying excitations and end up with an 
asymptotic form of the propagator which is valid for energies close to the ground state 
energy : 

G (t) ~ e ieBt det [ f; e ^f-kl)(t+i^)/2 X P X p' + ^-'^2)^^] (34) 

V=-oo / 

The threshold takes place at an energy u = 2e# — (2/7r) Jq f 5(e)de and the exponents 
are the same as for the repulsive potential. 

There is also a secondary threshold; we can find an asymptotic form of the propagator 
taking into account all the excited states in which the bound state is absent. The threshold 
takes place at an energy u = 2e F — (2/n) Jq f 5(e)de and the exponent is 1 — f3 + — f3~ with 
(3+ = (5 F /n - l) 2 and /T = {S F /n) 2 

V. CONCLUSIONS 



In this article we calculated the exponent of the spectral function in the asymptotic 
low energy range using the BA wavef unctions. The spectral function has no quasiparticle 
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peak and its incoherent part has a power law divergence at threshold. It is interesting 
to note that this singularity is accompanied by an orthogonality catastrophe with similar 
exponents; although such a coincidence was expected from perturbative analyses, no exact 
relation between these two quantities exist. These features are reminiscent of the static 
impurity problem although the exponents are different. This similarity however does not 
mean that the particle is localized in our case; in fact the effective mass is finite fo finite U 
]7]||. Care must be taken in distinguishing localization as probed by transport or optical 
experiments. 



Recently, Ogawa et al |L2] calculated the absorption spectrum for a Tomonaga-Luttinger 
model; they claimed that the exponents of the singularity did not depend on the mass of 
the particle; this was also valid for the spectral function. However our results show that 
the exponents differ in the two extreme cases = oo and = 1; in fact the authors 
neglected the backscattering and this process becomes relevant in the infinite mass limit as 
we already pointed out in section 3. We cannot exclude however that their results might be 
correct for any finite mass because the infinite mass problem is very particular due to the 
broken translation symmetry. A further analysis of the problem for intermediate masses is 
then needed. As no exact solution exist for other masses, we are extending our study of the 
spectral weight using Quantum Monte Carlo techniques. 



We mention the agreement of our results with the calculation in Ref. |L4] using conformal 
invariance : we recover their exponents in the presence of a strong magnetic field where the 
ground state of the full Hubbard model is almost ferromagnetic. Our results are also in good 



qualitative agreement with numerical results on this model [f[q| ; however no quantitative 
comparison can be performed because of the small systems studied numerically. 

Eventually, we can draw some conclusions concerning the experiments performed on 
quantum wires. Although we didn't calculate the absorption spectrum, we have noted the 
similarity of our situation with the core level problem and in the latter, the divergence 
of the spectral function is closely related to the edge singularities; thus the divergence of 
the spectral function in our model is consistent with the interpretation of the experimental 
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spectra in terms of excitonic effects. The calculation of the absorption spectrum with the 
BA solutions is under investigation. 
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VII. APPENDIX 

We prove here the assertions (11-13) which constitute central results of this paper. Sup- 
pose that we have found a set {kj} which are solutions of the B.A. equations. We first 
evaluate the overlap between two plane waves built up with different momenta k\ and k 2 . 

f L e m^)^s 2)dx = e ^- fa )^- fa)L - l = _ 4 sin(^)sin(^ 2 ) 
Jo i(ki — k 2 ) U 

where we have used the following property of the B.A. solutions : 

1 - e i{kl - k ^ L _ i(h - k 2 ) 
(1 - e ifclL )(l - e ~ ik2L ) ~ U 

We turn now to the evaluation of the normalization factor A, : 



(36) 



If 4sin(^) 2 _ smjSj) ~' (A + 1) sin(^) 2 / J_ 
M LU E/sin^) + (Eisin^,)) 2 U 2 



The sum Ezsin(^) is of order L and we recover the result (11). In order to gain insight in 
the scaling behaviour of this quantity, we can evaluate the sum for the ground state : 

l r kp . ( fU 



y^sin(5;)= [ sin farctan ( — 

i 7T Jo V \2k 



dk 



LU , / 2k F 
■ log -f + 



'2k F s 



+ 11 (38) 



We note that the corrections increase with U and inversely decrease with the density of 
fermions. This is generic for all the scaling behaviours we studied in this model. 
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We perform now the overlap between two functions <f>j and (j) n with j ^ n : 

= A,A n { ^ S f n ^ - L(SiD g^f (39) 

As the normalization factors behave like O (l/y/Z^j these overlaps decrease as well like 1/L. 

We end up the discussion with the asymptotic form of the Fourier components of our 
functions 4>j which are essential in the spectral analysis : 

— f\^*d,.(x)dx = AiVZf * n{Si) ( - 

y/L Jo 7T \(nj - to) - Sj/n 

1 E 7 S[n [ 6p) x , ) (4°) 



Ejsin(^) ^ (rip - to) - 5 p /ir, 

If to is close to the Fermi surface at uf = (N — l)/2, the last sum gives a logarithmic 
correction : 

j (nj - to) - di/7r 

The correction is then of order log iV/iV which decreases slowly and tends to shift the 
scaling region to larger sizes in any evaluation of the spectral properties for finite systems. 
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FIGURES 

FIG. 1. Exponent (3 + + (3~ of the orthogonality catastrophe as a function of the dimensionless 
parameter Uuf where U is the interaction strength and np the density of states at the Fermi 
energy for the equal masses case (solid line) and the infinite mass case (dashed line). 
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